Cyclotron resonance and Faraday rotation in graphite 
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The optical conductivity of graphite in quantizing magnetic fields is analytically evaluated for 
frequencies in the range of 10-300 meV, where the electron relaxation processes can be neglected 
and the low-energy excitations at the "Dirac lines" are more essential. The conductivity peaks 
are explained in terms of the electron transitions in graphite. Conductivity calculated per one 
graphite layer tends on average to the universal conductivity of graphene. The (semi)metal-insulator 
transformation is possible under doping in high magnetic fields. 

PACS numbers: 71.20.Di, 76.40.+b, 78.30.-j 



Graphite is usually considered as a layered semimetal 
composed of the graphene monolayers. Within this as- 
sumption, the graphite electron spectrum was evaluated 
many years ago within the Slonczewski-Weiss-McClure 
(SWM) theoryi. The so called "Dirac cone" turns into 
four bands with the twofold degenerate zero mode where 
electrons and holes are located^. The zero-mode disper- 
sion is very small, of the order of 20 meV, in the main-axis 
direction because the interlayer interaction is weak. Un- 
usual properties of graphite have attracted much atten- 
tion for more than 50 years. The most accurate method 
to study the band structure of graphite is a study of 
the Landau levels (LLs) through experiments such as 
magneto-optics^"^ and magnetotransporti^— . However, 
the interpretation of the experimental results involves 
a significant degree of uncertainty since, as noticed by 
Doezema et and Chuang at aS^ " it is not clear where 
the resonances are to be marked." 

The SWM theory requires the use of many tight- 
binding parameters and provides the simple description 
of observed phenomena either in the semiclassical limit 
of week magnetic fields or for high frequencies when the 
largest tight-binding parameter 71 = 0.4 eV plays the 
leading rolei2.. At the relatively strong magnetic fields 
B ^ 1-30 T and frequencies w ~ 10-350 meV, the smaller 
tight-binding parameters 72, 75, and A of the order of 20 
meV are essential. In this case, any physical property for 
graphite in magnetic field is represented by an integral 
over the momentum projection fc^. The SWM model can 
be simplified assuming that only the integration limits 
produce the main contribution o^i 1 5 , 1 7 _ g^Qj^ approxima- 
tion is similar to the theory of magneto-optical effects 
in topological insulatorsi^ and graphene^^. However, in 
the 3d systems, the other features such as the band ex- 
trema or the integration limits at the Fermi level can 
contribute as well. Therefore, the analytical expression 
for the dynamic conductivity in the presence of magnetic 
fields is needed for an interpretation of magneto-optics 
experiments. The theoretical study of magneto-optical 
properties in multilayer graphene is realized in Refi^*'. 

In this paper, we evaluate a formula for the optical con- 
ductivity of graphite in the presence of quantizing mag- 
netic fields and results are compared with experiments. 



We remind the notation for the LLs in graphite using 
the Hamiltonian in the form of Refs.^*''^^. The expres- 
sion for both the longitudinal and Hall dynamical con- 
ductivities is given. The (semi)metal-insulator transition 
is discussed in conclusions. 

Neglecting the trigonal warping 73 , the effective Hamil- 
tonian near the KH— line of the Brillouin zone can be 
written in the form 
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where k± — ^ik^ — ky, v — 1.02 x 10* cm/s is the intra- 
layer velocity, and jj are the functions of fc^: 

71 = 271 cos (fc2Co) , 72 = 272 COS (2fc2Co) , 

75 = 275 cos (2fc2Co) + A 

with the distance cq = 3.35 A between layers in graphite. 

At the magnetic field B, the momentum projections 
kx^y become the operators with the commutation rule 
fc„} = —2ehB/c, and we can use the relations 



k+ = y/2\e\hB/ca, fc_ = y/2\e\hB/ca+ , 

involving the creation and annihilation operators. We 
seek the eigenfunction of the Hamiltonian in the form 
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where the Landau number n > 2 and {x) are orthonor- 
mal Hermitian functions with n > 0. Then, every row in 
the Hamiltonian ([T]) becomes proportional to the Hermi- 
tian function with the definite n, and we obtain a system 
of the linear equations for the eigenvector Cs„ 
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FIG. 1: (Color online) LLs Ssn for n = to 6 in four bands 
s as functions of momentum kz along the KH— line in the 
Brillouin zone {K = 0, H = n/2co) at the magnetic field B 
= 20 T using the SWM model with the TB parameters taken 
from Refs.^i^ and neglecting trigonal warping; the Fermi 
energy ef = eV shown in dash-dotted line. Electron tran- 
sitions are possible between the levels £3,1 < and Es'n' > 
with the selection rule An = ±1 . 



where ujc = v 

At each Landau number n, the eigenvalues of the 
Haniiltonian are marked by the index of the band s — 
1, 2, 3, 4 numerating the levels from the bottom. We will 
use the notation \sn) for levels. LLs with n > 2 are in 
every four bands, shown in Fig. [l]as functions of along 
the KH— line of the Brillouin zone. The eigenvalues are 
determined by the equation 

[(75 -£)(72 -e) 1)] 
X [(75 - £)(72 - e) - ^In] ~ 71^(72 - e)^ = . 

In addition, there are four levels. One of them, Eq =72, 
with n = and the eigenvector Cq = (0,1,0,0). This 
level intersects the Fermi level and belongs to the electron 
(hole) band near the K (H) point. Other three levels, 
s = 1,2,3, are indicated hy n — 1 with C^i — 0. One 
level, |21), is very close to the level with n — 0. This 
pattern is consistent with Ref.— . The level structure at 
the K' H'—line is similar, therefore each level is fourfold 
degenerate, twice in spin and twice in pseudo-spin of the 
KH and K'H' valleys. 

In the region, where 71 ^ 72j 75i the two closest bands 
are written as 



£2,3 ('^) = 72 ± ujly/n{n - l)/7i . 

The Fermi energy at zero temperature, T—0, is de- 
termined by equality of electron and hole concentrations 
given by the sum of integrals 
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over the region in the Brillouin half-zone z = cok^ where 
electrons and holes are located. As known, the Fermi 
energy oscillates in the quantum Hall regime. We will 
consider the relatively strong magnetic fields where the 
oscillations do not exceed of 2 meV according to the 
electro-neutrality condition 

At finite temperatures the conductivity is expressed in 
terms of the correlation function^^ (for details see Ref.^^) 

V{uj)=tY^ I dxdx'Tr {v'-g {uj+,x,x')v^g {oj-,x' ,x)] , 

.(5) 

where Q (w-j., x, x') is the temperature Green's function, 
ui± = w ± LOm , Tr is taken over the |sn)— eigenstates, 
and the x, a;'— integration is over the coordinate involved 
explicitly in the Hamiltonian while the Landau gauge 
is used. Then, the Fourier transform of the eigenstates 
with respect the y— coordinate is assumed. Using the 
eigenfunctions ([2]), we write the Green's function of the 
Hamiltonian ([ij 



g'^P{u:,x,x') = Y^ 



The intralayer- velocity matrix is given by the deriva- 
tive of the Hamiltonian ([T]) 
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The straightforward calculation of the dynamical con- 
ductivity gives 
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where A, 



Y \Ci2 f^l* 1 f^i* r 

^ sn'^ s' n+1 ~^ ^ sn'^ s' n+l\ ' 

= Esn — £s',n+i IS the Icvel spacing, lOc 



v^2\e\hB/c is the cyclotron frequency, and /(e) = 
[exp(^^) -1-1]^^ is the Fermi-Dirac function. The in- 
tegration over the Brillouin half-zone, < z < 7r/2, and 
the summation over the Landau number n as well as the 
bands s, s' should be done in Eq. ([7]). 

The selection rule An = ±1 appears as a result of inte- 
gration over X and x' in Eq. ^ . If the trigonal warping 
is taken into account, the selection rule is change d^'^'^^ . 
Let us notice that the choice of the selection rule can 
be done examining the intensities of the cyclotron reso- 
nance lines. Our choice corresponds with observations in 
Refsiiii^. The results for the longitudinal conductivity 
ax'xiijj) are shown in Figs. [5] and 2] for two values of the 
magnetic field in order to compare the effect of the field. 
The dynamical Hall conductivity axy{uj) shown in Figs. 
13] and [5] describes the Faraday rotation^. The formula ([7]) 
is valid in the coUisionless limit, when the relaxation rate 
is much less than the frequency, F ^ a;. 
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FIG. 2: (Color online) Real and imaginary parts of the longi- 
tudinal dynamical conductivity at B =20 T (thick line); the 
partial contributions of various electron transition are shown 
in the thin lines. Temperature T=0.1 meV is less than the 
level broadening F = 5 meV. 
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FIG. 3: (Color online) Real and imaginary parts of the dy- 
namical Hall conductivity at B —20 T (thick line) with the 
partial contributions of various electron transition shown in 
the thin lines. 



In the calculations, we use the values 71 = 400 meV, 
72 = -10 meV,75 = 5 meV, A = 45 meV of Reis,^^, 
and r = 5 meV, Ep = S meV. 

The conductivity units of 



Ahco 



have the simple meaning, being the graphene dynamic 
conductivity e^/4?}^^ multiplied by the number I/cq of 
layers within the distance unit in the z— direction. One 
can see in Figs. [5] and HI that the value of conductivity 
calculated per one graphite layer tends on average to the 
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FIG. 4: (Color onhne) The same as in Figs. [2] but at B =7 
T. 
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FIG. 5: (Color online) Real and imaginary parts of the dy- 
namical Hall conductivity at B —7 T (thick line) with the 
partial contributions of various electron transition are shown 
in the thin lines. 



graphene universal conductance. 

Let us analyze the spectroscopy of the cyclotron res- 
onances at i? = 20 T. The line at 292 meV is a dou- 
blet resulted from the electron transitions |25) — >■ |36) 
and 1 26) — )• |35) near the K point of the Brillouin zone 
(see Fig. [IJ. In the similar way, the 240- meV line in- 
cludes the 1 24) 1 35) and |25) |34) doublet sphtted 
due to the electron-hole asymmetry. However, the broad 
line at 165 meV involves other transitions besides the 
similar |23) ^ |34) , |24) ^ |33) doublet at 184 meV. 
First is the transition |11) — > |10) (161 meV) near the 
H point. Then, the transitions |22) -> |21) produce the 
broad band. The high-frequency side of the band (170 
mev at H) and the low-frequency side (70 meV, at the 
intersection of the |21)— level with the Fermi level) con- 
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tribute into the 165-meV and 50-meV lines, correspond- 
ingly. The position of the broad band is very sensitive to 
the variation of the tight-binding parameters and to the 
magnetic field. The comparison of Fig. |3]to Fig. [2]shows 
that the |23) |34) line does not interfere with other 
transitions at the weaker field. The main contribution 
into the sharp 50-meV line is resulted from |21) — )■ |32) 
transitions near the /("-point. The positions of the lines 
for the fields in the range of 20 - 30 T agree very well 
with observations of Refsjiii^. We do not correct the posi- 
tions by a variation of the tight-binding parameters. Let 
us emphasize that the imaginary part of the dynamical 
conductivity is of the order of the real part. 

The optical Hall conductivity axy{uj) in the ac regime 
is shown in Figs. [3] and [H It is evident that the inter- 
pretation of the Faraday rotation governed by the con- 
ductivity <Txy{oj) (see Fig. ??) is much more complicated 
in comparison with the longitudinal conductivity. The 
conductivities axx{^) and axy{^jj) allow us to calculate 
the reflectivity and the Faraday rotation as functions of 



frequency. 

Notice that the band structure shown in Fig. [1] 
constrains to consider the (semi)metal-insulator transi- 
tion while varying the carrier concentration and apply- 
ing the magnetic field. The phase transition induced 
by the electron-hole interaction has been discussed in 
Refs i^°i^^'^^ . As one can see in Fig. [T] the hole dop- 
ing can decrease the Fermi energy. While the Fermi level 
appears between the |22) and |10) levels, the insulator 
arises with a gap of 34 meV. This phase transition does 
not involve the electron-electron interaction and it is re- 
sulted from the layered graphite structure with the small 
electron dispersion of the zero mode in the fc^-direction. 

Author is thankful to A. Kuzmenko for discussions and 
corrections. This work was supported by the Russian 
Foundation for Basic Research (grant No. 10-02-00193- 
a) and by the SCOPES grant IZ73Z0.128026 of the Swiss 
NSF. Author is grateful to the Max Planck Institute for 
the Physics of Complex Systems for hospitality in Dres- 
den. 



^ J.W. McClure, Phys. Rev. 108, 612 (1957); J.C. Slon- 

czewski and P.R. Weiss, Phys. Rev. 109, 272 (1958). 
^ G.W. Semenoff, Phys. Rev. Lett. 53, 2449 (1984). 
^ W.W. Toy, M.S. Dresselhaus, G. Dresselhaus, Phys. Rev. 

B 15, 4077 (1977). 
■* R.E. Doezema, W.R. Datars, H. Schaber, A. Van Schyndel, 

Phys. Rev. B 19, 4224 (1979). 
^ Z.Q. Li, S.-W. Tsai, W.J. Padilla, S.V. Dordevic, K.S. 

Burch, Y.J. Wang, D.N. Basov, Phys. Rev. B 74, 195404 

(2006). 

^ M. Orlita, C. Faugeras, G. Martinez, D.K. Maude, M.L. 
Sadowski, M. Potemski, Phys. Rev. Lett. 100, 136403 
(2008). 

'' M. Orlita, C. Faugeras, J.M. Schneider, G. Martinez, D.K. 

Maude, M. Potemski, Phys. Rev. Lett. 102, 166401 (2009). 
* M. Orlita, M. Potemski, Semicond. Sci. Technol. 25, 

063001 (2010). 

® I. Crassee, J. Levallois, A. L. Walter, M. Ostler, A. Bost- 
wick, E. Rotenberg, T. Seyler, D. van der Marel, A. Kuz- 
menko, Nature Physics 9, 48 (2010). 

1° Y. Kopelevich, J.H.S. Torres, R.R. da Silva, F. Mrowka, H. 
Kempa, P. Esquinazi, Phys. Rev. Lett. 90, 156402 (2003). 

" I. A. Luk'yanchuk, Y. Kopelevich, Phys. Rev. Lett. 97, 
256801 (2006). 

Z. Jiang, Y. Zhang, H.L. Stormer, P. Kim, Phys. Rev. Lett. 

99, 106802 (2007). 
" J.M. Schneider, M. Orlita, M. Potemski, D.K. Maude, 

Phys. Rev. Lett. 102, 166403 (2009)^ 

" A.N. Ramanayaka, R. G. Mani. ,.arXiv:1010.0603k l (2010). 



K.-C. Chuang, A.M.R. Baker, R.J. Nicholas, Phys. Rev. B 

80, 161410(R) (2009). 
" L.A. Falkovsky Phys. Rev. B 82, 073103 (2010). 
^"^ G. Li, E.Y. Andrei, Nature Phys. 3, 623 (2007). 

W.-K. Tse, A.H. MacDonald, Phys. Rev. B 82, 161104R 

(2010). 

T. Morimoto, Y. Hatsugai, H. Aoki, Phys. Rev. Lett. 103, 
116803 (2009). 

2° M. Koshino, T. Ando, Phys. Rev. B 77, 115313 (2008). 
A. A. Abrikosov, L.P. Gorkov, I.E. Dzyaloshinski, Methods 
of quantum field theory tn statistical physics, Dover Publi- 
cations, N.Y. 

L.A. Falkovsky, A.A. Varlamov, Eur. Phys. J. B 56, 281 

(2007) . 

P. Nozieres, Phys. Rev. 109, 1510 (1958). 

A. B. Kuzmenko, E. van Heumen, F. Carbone, D. van der 
Marel, Phys. Rev. Lett. 100, 117401 (2008). 

H. Suematsu, S. Tanuma, J. Phys. Soc. of Jpn. 33, 1619 
(1972). 

B. Partoens and P.M. Peeters, Phys. Rev. B 74, 075404 
(2006). 

A. Griineis, C. Attaccalite, L. Wirtz, H. Shiozawa, R. 
Saito, T. Pichler, A. Rubio, Phys. Rev. B 78, 205425 

(2008) . 

D.V. Khveshchenko, Phys. Rev. Lett. 87, 246802 (2001). 
2^ E.V. Gorbar, V.P Gusynin, V.A. Miransky, L.A. 
Shovkovy, Phys. Rev. B 66, 045108 (2002). 



